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Abstract 

A general approach for evaluating the nonfactorizable contributions to the 
hadronic matrix elements is described in terms of the higher twist wave func¬ 
tions of the mesons. An example is illustrated for the B° —> 7r + 7r“ decays. 
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The factorization approximation has conventionaly been used to describe 
the two-body nonleptonic decays of D and B mesons. However, it has been 
observed [jlj that the naive factorization approach to nonleptonic heavy me¬ 
son decays is incompatible with experimental data. The available data for 
D-meson decays seems to support the so-called rule of discarding 1/N C cor¬ 
rections. This rule has been formulated in ref. ||, which is in contrast to the 
standard prescription that keeps the l/iV c term. The dynamical explanation 
of this approximate cancellation of 1/N C terms in the D decays has been pre¬ 
sented in ref. ||. From such an explanation, it appears that discarding l/iV c 
terms may not necessarily be a universal rule. The recent experimental data 
on B-meson decays suggests that discarding l/iV c terms is not a universal 
rule for all channels. To understand these problems, one needs to see what 
approximation is being made in the naive factorization approach and how 
to systematically improve upon the approach. It will be seen that the issues 
concern the nonfactorizable contributions of the operator of two color-octet 
currents. 

The leading nonfactorizable contributions to hadronic matrix element 
arise from the soft gluon interactions. It turns out, such soft gluon con¬ 
tributions to hadronic weak decays of mesons can be further factorized in 
terms of twist wave functions (TWF) of the mesons. We shall describe in 
this note a general approach for their evaluation. 

As an illustration, let us consider an effective QCD-corrected weak Hamil¬ 
tonian with two familiar operators 0\ and C> 2 - 

Heff ~ Ci(/l)Oi + C2(/x)02 ( 1 ) 

where Cj(/r) are the Wilson coefficient functions and O* the effective four 
quark operators. For £> —> 7T7 t decays, Oi are given by 

O i = (ub)v- A (du) V - A , 

C >2 = ( U a bp)v-A(di3U a )v-A (2) 

In the factorization approximation, the hadronic matrix element is factor¬ 
ized into the product of two matrix elements of singlet currents, which are 
evaluated in terms of decay constants and formfactors, and characterized by 
the parameters and These parameters are related to the Wilson 
coefficients c* 
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( 3 ) 


(o) i (o) ■ 

Ol =Ci + —c 2 , «2 =C 2 + — Cl. 

This naive factorization fails to describe the nonleptonic charm meson 
decays. The approximation neglects the contributions from the product of 
two color-octet currents, 

Of = 2 (dt a b)v-A(ut a u)v-A, Of = ( ut a b)v-A(dt a u)v-A (4) 

Where t a (a — 1, • • •, 8 ) are SU(3) C generators, and the color indices in the 
operators Of have be omitted. It becomes manifest that in the naive factor¬ 
ization approach at this level the contribution to the B —> n + TT~ decay from 
the operators Of vanishes due to the mismatching of color. 

Using the algebraic relation = ^(5 aa '8gg> — jr8 a p5 a >p>), and ap¬ 

plying the Fiertz transformation, one easily obtains the following identities 

o, = 2-o, + Of (5) 

where (i,j) represents the pair (1,2) and N c is the number of colors. The 
leading nonzero contribution of the operator Of to the B —» decay 

arises from the soft gluon emission which recovers the color matching of the 
octet currents to the hadron. The resulting new effective operators due to 
the soft gluon interaction can be written as a product of two color-singlet 
currents again. The leading nonzero contribution is given by 

A 0 c — —i J dx < n + (p_)\T (j(x)Of(0))\B(p) > (6) 

where j(x) = d[x)A(x)d{x) with A{x) = g s A{x)f l t a ^^’ (other quark flavors 
contribute to j(x), and their effects can be similarly calculated). Contracting 
the down quark pair, and using the algebraic relation of the SU(3) generators 
t a t b = jf8 a b + |(/ a6c + id abc )t c , one then obtains 

1 c r dn p iQx 

M = Y C J dx J ( 4 7t)4 q * - ml < > (7) 

with 

J^ = d(x)Y(q P l P + m d ) 7^(1 - 7 5 )«(0), 

J„ij,= [ da a x p u(0)G pu (ax)'y fl (l - 75 ) 6 ( 0 ). ( 8 ) 

J 0 
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where we have neglected the terms given by two color-octet currents since 
their nonzero contributions involve the next-to-leading order contributions 
of the coupling constant g s . 

At this level, one can again use the factorization approach to evaluate the 
hadronic matrix element, and obtain 




As the currents q) and J^ v are nonlocal, their hadronic matrix elements 

will be evaluated by so-called twist wave functions 



•[(jA + V-Q v - Q ■ P-9^ - ie VW(J qpP-a) 0 a («) (10) 



where <f>A , 4>p and are pion twist wave functions defined as follows^, H 




( 11 ) 


< 7T (p_)|M(x)cr All/ 7 5 d(0)|0 > 


( 12 ) 



(13) 


The twist wave functions (j)p, (j>T and 4>a satisfy the normalization conditions 



(14) 


These wave functions can be obtained from QCD sum rules. The results 
depend on the renormalization group scale. For p 2 = m 2 B — m 2 = 5 GeV 2 , 
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they are given by[Q: <j>A(ot, /i) = 19.32a(l — a)[(2a — l) 2 + 0.11], 4>p(a, /i) = 
f[1.078 + 0.766(2a- l) 2 ], 0r(a,/x) = a(l - a). 

To find the matrix element < 7r + (p + )| J^ \B(p) >, let us consider the 
correlator 


1 ?,m(P+,P,x) 


i J dye tpy < tt + (p + )\J^ b(y)i^ 5 d(y)\0 > 


m%f B 1 
m b m B — p 2 




(15) 


where < B(p)\b(y)i'y 5 d(y)\0 >= (m 2 B f B /m b )e^ ipy . The ellipses denote the 
excited and continuum states which are usually taken into account by the 
dispersion integral with the effective threshold. The single soft gluon contri¬ 
bution is given by 

r r e ~i(p~k)y /-l 

T=(P + ,P,I) = i J d Vj dk (4*)*(ml-k*) Jo da ° a ° X ‘’ (16) 

• < n + {p + )\u(())Gp U (aox)^^(k a ^ a + m b )^d(y)\$ > 


which can be evaluated by introducing the pion wave functions with higher 
twist 

T ">+' P ' *) =>Jd„j dk (i7r) K _ t!) J o 

US 3 ,J DufeWe**'”’**”*! 

■[{k v x ■ p+ - p +u x ■ k)p + fj - (. g Vfl x ■ p + - p+uX^k ■ p + \ (17) 

+fAm b J Da i e w+{o ‘ 2 y+a 3 ax) [(f) ± (a i )(g uti x ■ p + - p+^x^) 

,, , . . , ... y ■ x x ■ p + w , m 2 

+ ( 0 ||(«i) + (j)±{a i )){p +u - VvAP+p - a 2 k„ -) j 

y ■ P+ y ~ P+ m d rn b 

where </> 37r and (p\\, (ft± are the pion wave functions of twist 3 and 4, respec¬ 
tively. They are defined as || U 

< 7r + (g)|h(0)Gp^ (T 75(i(?/)|0 >= -i J Dai(j) 3 n (ai)e iq ( a2y+Q3Q ° x) 

■ fn^ 3 ir\Qfi{dp 9 iycr Qvdpcr) da^QpQvp Qv 9 pp)\i (l§) 

<ir + (q)\u( 0 )Gp^p, 7 5 d(y)\ 0 >= J Da^a^ 2 ^ 3 ^ 
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'fir^4n ' \dlv{.9pfi 


Vp% 

y-Q 


) Qpidvfj, 


Vv% 

y ■ q 


+/X^ Wl- W / Ba^||(a,:)e‘’ ( “ a ' + “»“" I) 


( 19 ) 


where Dai = da\da 2 da^8(a\ + «2 + «3 — 1). The constants <5 3?r and 5^ are 
renormalization scale dependent and have the values|5|] <5 37r ~ 0.027GeV and 
5 47r ~ 0.45 GeV at n ~ 1 GeV. The results for </> 37r , (fi± and 4>\\ are given 
by </ , 37r( Q; i) = 360aia 2 c> ! 3 [l T ^i,of (7cr 3 — 3) + c<^2,o(2 — 4oqa 2 — 8<a 3 + 8a 3 ) + 
uji t i[?>aia 2 — 2cr 3 + 30(3) + •••], 4>±.(oti) = 30(aq — cx2)cr 3 [^ + 2£'(1 — 2<a 3 )], 
4>\\ (cq) = 120e(ai — a 2 )aia 2 a 3 • where a^o = —2.88, a; 2) o = 10.5, = 0 

and e ~ 0.5 at p ~ 1 GeV. 

The integral of eq. (18) for the variables y and k can be performed and 
is given by 


T vll (p+,p,x) 


= i Daj 


u 


da 0 a 0 e iaoa3P+a 


m 2 b - (p - a 2 p + ) 2 ) Jo 

{537T037T {0ti){p v X ■ P+ - P+yX ■ P)p +ll - {g vlJ X ■ P + - P+yXfj) 

■ ({p ■ p + - a 2 p+)^0 37r (ai) - Sl w m b ^±(ai)^ (20) 

+Slmb ( Mai) + *> (ai) ) r^ m |’-(p-(a 2 T/3K) 2 

ml. 2 

•- (9vpX ■ p + -p +v x li ) + —- - - - - ——— 

m d m b mi - {p - ( a 2 + p)p+) 

\p u x ■ p + - p +u x ■ p) (p +li + -^^a 2 ((a 2 + (J)p+p +Pp) ]]} 

V m d m b I 


Here we have used the relation 1 /p + ■ y = —i / 0 °° d/3e l/3p+ ' y . With the above 
results, we further consider 
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dq e iqx 
(47t) 4 q 2 — ml 


T v ^{p-,q,x)T^{p + ,p,x) 


m B f B 1 A , f°°p(ml,s)ds 

- 2 - 2 A °2 + / -o 

m b m B — p z 2 Js 0 s — p z 


( 21 ) 


where the dispersion integral with the effective threshold So assembles the 
excited and continuum states. By substituding eqs.(10) and (24) to the left 
hand side of the above equation, it is not difficult, though tedious, to integrate 
out the variables x and q. The result is found to be 
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1 

iVc 

= i- 


dx 


dq 


rj—iqx 


N r 


(47r) 4 q 2 — m 2 d 

fl [ da f da 0 a 0 f Dai —— 

Jo Jo J mi — [p — i 


T u ^{p-,q,x)T^(p + ,p,x) 

{ 


P 


l-ip- «2 p+y 

a(a 0 a 3 p 2 + 2(a - a 0 a 3 )ml ) 


2 [(a - a 0 a 3 ) 2 7nl + aa 0 a 3 p 2 - m 2 ^ 2 
■{d 3n (f) 3v {ai)p 2 - 5l n m b (j)±(ai) + + (p\\ (a*)) 


dfd 


m 2 b - {p- a 2 p+y 


p 


o m 2 b - {p- (a 2 + f3)p+) 2 ml-{p- (a 2 + /?)p+) 2 


( 22 ) 


m; 


- 0^2 (1 + Ol2 + /?)) — 


m; 


-)) 


'(1 + 

m d m 6 

+ ((a - a 0 a 3 ) 2 m 2 + aa 0 a 3 p 2 ) (<5 37r 0 37r (ai)p 2 - 2<f 2 7r m b 0_ L (a i ) 


- 254 7 r m b ( 0 _ L (aj) + </>||(cq)) 

p 2 


d/3- 


m b ~iP~ &2P+f 


mt 


)] 


(a - a Q a 3 ) 2 m 2 + aa 0 « 3 P 2 - m d 

/»oo 

+5f jr m b (0j.(a i ) + 0||(a*)) / 


o ^ ml - ip - (a 2 + /?)p+) 2 m d m 6 
[25 37r 0 37r (aj)p 2 - 354 7r m 6 0 ± (a i ) 


ml - ip- a 2 p+f 


o ml-ip - (a 2 + /?)p+) 2 


(1 + 


m 


m d m b 


-a 2 {l + a 2 /3)) 


p 


mz 


ml - ip - (a 2 + /?)p+) 2 m d m b 


]} 


where we have neglected the terms of order m 2 /m 2 B . The remaining integral 
may be carried out numerically (in fact, the integral for the variable (3 can 
be explicitly integrated out). Before doing so, one may use the QCD sum 
rule approach^ by applying the Borel transformation in the variable p 2 , i.e., 


L M ifiQ 2 ) 


lim 

Q 2 ,n—>oo,(Q 2 /n)=M 2 


(g 2 )(n+i) 


d 

dQ 2 


fiQ 2 ) = /(M 2 ) 


(23) 


to the left and right hand side of eq.(25) with use of the more explicit ex¬ 
pression eq.(26) for the left hand side. M is the so-called Borel parameter. 
It can be seen that the Borel operator transforms the series in Q~ 2 into a 
series in M~ 2 : 


6 

















n 


(24) 


From duality considerations, the dispersion integral on the right hand 
side of eq.(25) is cancelled against the corresponding part of the dispersion 
integral of the QCD result on the left hand side. One then reads off the 
needed amplitude A 0 c. The procedure is quite standard following the usual 
QCD sum rule calculationslJ. The utility of the result is that it depends on 
twist wave functions that can be calculated from two point Greens’ functions. 
We have used B° —> it + tt~ as an example but the method applies equally 
well to other heavy to light channels. A detailed application will be discussed 
elsewhere. 

In conclusion, we have shown that the nonfactorizable contributions to 
the hadronic matrix elements of the nonleptonic two body decays can be 
calculated directly from the QCD theory. It has been seen that the leading 
nonzero contributions from the single soft gluon emission can actually be fac¬ 
torized again in terms of the higher twist wave functions of the mesons. It is 
expected that the approach presented here will be helpful in the understand¬ 
ing hadronic decays of the heavy mesons, especially for the heavy to light 
meson decays. It shall provide a complement to the approaches presented in 
ref. (Tj for the calculations of the B —> Dn and B —> J/i^K decays in which 
one of the final hadrons is heavy. 
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